Dynamics of genuine multipartite correlations in open quantum systems 
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We propose a measure for genuine multipartite correlations suited for the study of dynamics in 
open quantum systems. This measure is contextual in the sense that it depends on how information 
is read from the environment. It is used to study an interacting collective system of atoms undergoing 
phase transitions as external parameters are varied. We show that the steady state of the system 
can have a significant degree of genuine multipartite quantum and classical correlations, and that 
the proposed measure can serve as a witness of critical behavior in quantum systems. 
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Introduction. — There is a lot of interest in the char- 
acterization and quantification of correlations between 
quantum systems. The fact that quantum systems can 
be correlated in ways that classical systems can not is 
generally understood to be essential to emerging quan- 
tum technologies [l|, and thus understanding the na- 
ture of such correlations is one of the most active ar- 
eas of research in quantum information theory. Recently, 
there has been much effort made towards quantifying 
correlations — quantum or classical — for multipartite (as 
opposed to bipartite) systems. In particular, there exists 
a question of whether or not a system contains genuine 
multipartite correlations, i.e., are there correlations that 
cannot be accounted for by any subsystem alone. A gen- 
eral framework, where a measure for a type of correlation 
was defined as the "distance", as measured by relative 
entropy, from one state with a defining property to the 
closest state without said property, was recently proposed 
by Modi et al. in (2j. This framework permits a unified 
view of the measures of relative entropy of entanglement 
3, quantum discord Q, classical correlations @, quan- 
tum mutual information and several new quantities. A 
powerful feature of this framework is that the suggested 
measures are well defined, and relations between them 
also hold for multipartite systems, thereby allowing for 
natural extensions of the above quantities beyond the bi- 
partite case. A natural follow up to the work in Q was 
done by Giorgi et al. in [f| where n-partite quantum 
and classical correlations were further divided into gen- 
uine and non-genuine parts. This is not the only effort in 
this direction however, and it is not known how to cor- 
rectly measure genuine multipartite correlations. Inves- 
tigations of the validity of the various measures involving 
the dynamical evolution of correlations in experimentally 
relevant situations therefore seem beneficial. This is the 
purpose of the work presented here. 
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A notorious feature of many of the proposed measures 
of correlations in quantum systems is that they involve 
difficult extremizations over operators or states. This 
makes calculations hard in practical applications and re- 
stricts analytical results to simple examples. Even nu- 
merically, the extremization required is often intractable. 
Quite generically though, the extremization is not neces- 
sary for pure states, for which the measures are greatly 
simplified. The best known example is for bipartite en- 
tanglement, for which all the various mixed state entan- 
glement measures reduce to the entropy of the reduced 
state of one of the subsystems. This gives hope for a 
dynamical characterization of correlations when the dy- 
namics themselves naturally give rise to a decomposition 
of the quantum state into an ensemble of pure states. 
For bipartite entanglement this has been exploited to de- 
fine an average entanglement for open Markovian sys- 
tems with an evolution described by a Lindblad master 
equation Q, and much work has been done in recent 
years to investigate this entanglement measure in various 
experimentally relevant situations |8l-[l3|. We present an 
analogous definition for genuine multipartite correlations 
and show how it can be used to study the dynamical evo- 
lution of correlations in open multipartite quantum sys- 
tems. We explore this new measure by applying it first 
to very simple three-qubit systems, before moving on to 
the more physically interesting case of a collective spin 
state of a large number of atoms that undergoes phase 
transitions as an external parameter is varied. For this 
last case, the proposed measure of genuine multipartite 
correlations is shown to be a good witness for the critical 
behavior of the system. 



Average genuine multipartite correlations. — We use 
the definitions of genuine n-partite correlations put forth 
in 0. As was shown, genuine correlations can be con- 
sidered to be bipartite measures for some particular bi- 
partite cut of the system, and for a pure state p gen- 
uine classical and quantum correlations therefore coin- 
cide, and are given simply as the entropy of the reduced 
state of subsystem i, minimized over all possible biparti- 
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tions HQ: 

jW( p ) =5 W W = r W( p) /2 = m i n [ % )]. (i) 

z 

Here, t7^"\ 2?W and 7~' n ' refer to the genuine classi- 
cal, quantum, and total correlations, respectively. We 
note that this quantity was also considered as a multipar- 
tite entanglement measure in fl3| . We wish to consider 
the dynamical evolution of these correlations for an open 
quantum system consisting of several sub-systems. We 
will consider an environment such that the Born-Markov 
approximation applies and that the reduced density ma- 
trix of the system, p, evolves according to a master equa- 
tion in Lindblad form [14] (h = 1) 

p=-i[H,p]+Y J D[Ji]p- (2) 

i 

Here D[0]p = OpC - \O^Op - \pO^O, and the 
Ji are "jump operators" acting on the system, de- 
scribing the environment in terms of various types 
of noise. Given a pure initial state p = the 
quantum trajectory method [HI, [l6[ evolves the state in 
time step dt from \ip(t)) to either Ji\ijj(t)} /\\ Ji\ip(t)\\, 
with probability pi = (ip(t)\ jj Ji\ifj(t)) dt , or to 
exp(-iH eS t/h)\ii(t))/\\ exp(-iH eS t/h)\ij(t)) 1 1 , with 
probability 1 - Y^iPh where H cS = H - jjji- 
This means that the state \ip(t)) evolved in this way is 
pure for all times t, and Eq.([T]) applies. The density 
matrix p(t) can be retrieved as p(t) = \ipj(t))(%ljj(t)\, 
where each \ipj(t)) is a stochastic realization according 
to the quantum trajectory method, and the bar denotes 
an average over all possible realizations. We thus define 
the average genuine n-partite correlations for an initially 
pure n-partite state with time evolution given by Eq.Q 
as 

4 n) (<) = miii) - nun S(tv z \^(t))(^(t)\), 

(3) 

where as before i runs over all bipartitions, and tr, de- 
notes the corresponding partial trace. 

For bipartite pure states, entanglement represents all 
quantum correlations Q, and indeed we see that Eq.Q 
reduces to the average entanglement over trajectories de- 
fined in [tJ for n = 2. The symbol U here denotes the 
fact that this quantity is dependent on the choice of un- 
raveling of Eq.Q, i.e., the master equation is invariant 
under the unitary and inhomogeneous transformations 
^ G. t = J2j u rjJj + on, H ->■ H + l/2i(a*Ji - ctij}), 
where iiy is a unitary matrix and the Ui are complex 
numbers. In other words, there is a freedom in the choice 
of jump operators, and different choices will in general 
produce different trajectories and therefore a different 
time evolution for the average genuine correlations. 

Three-qubit systems. — To illustrate some of the new 
behavior in going from the bipartite to the multipar- 
tite case, we first concentrate on the simplest example 



of three qubits. We restrict our attention from the vast 
variety of different interactions between the qubits and 
their environments to the case of spontaneous emission 
of atomic qubits with the qubits self-Hamiltonian set to 
zero, H = 0. That is, we consider unitary combinations 
of the jump operators ^/tao, ^/tb^ and sflc c i where 
1a.b,c are decay rates and a = |0)(1|.a is the lowering 
operator on qubit A and similarly for b and c. 

As a first example, we consider the initial state = 
(2|011) + 2|101) + 1|110))/V9 and the case where just 
qubit B is decaying (i.e., jb > 0, ja — 7c = 0). 
The initial qubit entropies are S(pa) — S(ps) = 1 and 
C^ 3) = S(pc) = 1/2. Due to the decay of qubit B 
the entropy S(pB(t)) will eventually become smaller than 
S(pc(t)), so that Cjj ' has a smooth behavior everywhere 
except at the point in time when the entropy of qubit B 
is equal to that of qubit C . This is illustrated in the left 
panel of FigfT] This non-smooth behavior is only possi- 
ble for three or more qubits and ultimately stems from 
the symmetrization over bipartite cuts in the definition 
of genuine n-partite correlations. 

Next we consider an example where the coupling of 
the system to an environment can create genuine corre- 
lations even if no such correlations are present initially. 
We take \i>) 2 = 1/V2(|010) + |001>) as our initial state. 
This gives the initial entropies S(pb) — S(pc) = 1 an d 
S(pa) — 0. Under a direct photo-detection scheme, i.e., 
jump operators ^fjXa, ^fB~b and ^/tcc, the genuine 3- 

partite correlations C\j will stay zero for all times. If 
two of the decay channels are combined (e.g., in a beam 
splitter) however, we have the possibility of genuine 3- 
partite correlations being created in the system from 
measurement back-action. We consider the unraveling 
Ji = l/V2(y/7Xa + i^b), J 2 = l/V2(i y /y A ~a + ^yEb) 
and J3 = y / 7cc. In the right panel of FigfT] we show 

(3) 

Cjj for this unraveling as a function of 7^ for three 
different 73 (and 7c = J a), showing that the maximum 
of the genuine 3-partitc correlations in fact grows with 
the damping of qubit B. Similar counter-intuitive be- 
havior of environment-assisted entanglement that grows 
with the damping of a cavity mode has been shown previ- 
ously for bipartite entanglement 0] , but it has a different 
source in our case, namely that the entropy of a subsys- 
tem increases under a detection scheme where its output 
is mixed with another subsystem of higher entropy. 

Critical behavior in collective spin systems. — 
Collective spin systems are an interesting class of 
many-qubit systems for which states and dynamics 
are symmetric with regards to permutations of the 
qubits and can be described, e.g., in terms of Dicke 
states [l7| . which are simultaneous eigenstates of the 
collective angular momentum J and its z component 
J z . Physically, such systems arise when the dynamics is 
such that excitations cannot be uniquely identified with 
individual qubits. An A^-qubit symmetric Dicke state 
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FIG. 1: (Color online). Left: Average entropies for the three 
qubits, S(pa) (blue), S(pb) (green) and S(pc) (blue) for ini- 
tial state \ip)i and with decay of qubit B only. The minimum, 

(3) (3) 

C v , is marked by the thick black line. Right: Cjj for the 
initial state \ip)a when the decay channels of qubits A and 
B are mixed at a beam splitter before detection, for three 
different 7_g. 



with m excitations is given by 18 1 



\m,N) 



N 



-1/2 



^ Pk ( 1 1 1 , ■ • ■ , lm, m+ l, . . . , Oat)) , 

(4) 



k 



where {Pk} denotes the set of all distinct permutations 
of the spins. 

Quantum phase transitions have been studied theo- 
retically for such systems and shown to be marked by 
critical behavior of bipartite entanglement measures (be- 
tween either pairs of spins or blocks of spins) [18l422| . 
and recently in the quantum correlations between pairs 
of spins [23]. We show now that this is also the case 
for our measure of genuine iV-partite correlations. The 
specific (and topical) model we consider is the Lipkin- 
Meshkov-Glick (LMG) model Q, which is described by 
the Hamiltonian 



LMG 



-2hJ z -2(\/N){J 2 x+1 J 2 v ) 1 



(5) 



where the J x , J y and J z are collective angular momentum 
operators, h and A are effective magnetic field and spin- 
spin interactions strengths, respectively, and 76 [— 1, 1] 
is an anisotropy parameter, which we will set to zero. A 
dissipative version of this model was studied in (2l| based 
on the collective interaction of an ensemble of atoms with 
laser fields and field modes of a high-finesse optical res- 
onator. Under appropriate conditions, the cavity modes 
may be adiabatically eliminated to give a master equa- 
tion describing the collective atomic system alone, 



P = ~ «[#LMG, p] + JjD[J-{ 



(0) 



where J+ = J2j i s tne collective raising operator 

and J- its Hermitian conjugate. The quantum jumps 



described by J + can be associated with photo-detections 
of one of the cavity modes. 

Eq.© has a steady state solution, p ss , that exhibits 
both 1st and 2nd order phase transitions as the effective 
magnetic field strength h is varied in relation to the spin- 
spin interaction strength A; in particular, around h = 
and h = A(> 0), respectively [21j. To find the average 
genuine TV-partite correlations, Eq.©, we note that a 
pure state of the system, evolved according to the quan- 
tum trajectory method, can be written as a linear combi- 
nation of states of the form Eq.(U). Any such symmetric 
Dicke state has a Schmidt decomposition under the ar- 
bitrary bipartition (1,2,..., N^iNi + 1, Ni + 2, . . . , N) 
that can be written as [lg 



\m,N) = '%2\%\k,N 1 )®\m-k,N-N 1 ), (7) 

k 

where A™ = O'^)^ i^)^ • F <* a ^te 
P = Em n c m c n\ m i N}{n, -^1 we can then find the state 
reduced over the first Ni qubits: 

PN X v /' = c ™< J2 Ki^lT 

m,n k 

x \m- k,N - Ni)(n- k,N - A^|. (8) 

The genuine A-partite correlation is then the entropy 
of this state minimized over all possible bipartitions (p 
pure) : 



jW(p) = mm S(p Nl ). 



(9) 



In Fig(2]we show numerical results for the average gen- 
uine ./V-partite correlations, Eq.Q, after letting the sys- 
tem reach its steady state and averaging over a large 
number of quantum trajectories with the jump operator 
y / 2Tb/NJ + . The parameters used were Tb — 0.2 and 
N = 10, 25 and 100 qubits. With increasing N, the be- 
havior of points clearly to 1st and 2nd order phase 
transitions around h/X = and h/X = 1, respectively. 
In the inset of the figure we also show the average of 
maxATj S(pn 1 ). Although not a measure of genuine cor- 
relations, this quantity also shows very marked behavior 
around the critical points. 

We can also study alternative, diffusive unravelings 
of the master equation Eq.([S]) by considering, e.g., ho- 
modyne detection, in which the output cavity field is 
combined at a beam-splitter with a strong laser field 
before detection. In the present context, this corre- 
sponds to measurement of the spin quadrature Xg — 
c l6 J + + c~ lB J_. The evolution of a quantum trajectory 
\ipj(t)) can then be given as a continuous stochastic dif- 
ferential equation jlj, [25[ , 



d|^(*)) 
2N 



-iH\4>j{t))&t 



(X 9 



J-J+-\{Xgf }\, ;,(/))<!/ 




\^(t))dW(t), 



(10) 
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FIG. 2: (Color online) Average genuine iV-partite correla- 
tions, , versus h/X for the dissipative LMG model with 
7 = 0, r 6 /A = 0.2, and N = 10 (red), N = 25 (green) and 
N = 100 (black) qubits. In the inset we show the complemen- 
tary quantity that maximizes the entropy over all bipartitions. 



where dW(t) is a real Wiener increment. Results for 
C\f ^ with this unraveling are shown in Fig|3]for N = 100 
and two different detection angles 0. Considerable sen- 
sitivity of Cjf to 6 is apparent in the region between 
h/X = and 1. In this region, the system state during 
a quantum trajectory can approach a superposition of 
J x eigenstates of opposite signs (note that the spin Q- 
function computed from p ss has a two-peaked structure 
along the J^-axis [Hj]), provided the homodyne measure- 
ment is unable to resolve these eigenstates, i.e., provided 
the phase 9 ~ 90°, corresponding to measurement of J y . 
Such a superposition of J x eigenstates is strongly corre- 
lated, giving rise to a large value of c\f \ while a single 
J x eigenstate, towards which the state is projected by a 
measurement with 6 — 0°, has no genuine correlations 
and hence a negligible value of c\f ^ over much of the re- 
gion between h/X = and 1. We note that the detection 
angle that gives the maximum value of varies with 
both h/X and N, and that the sensitivity increases with 
N. 

We can of course write the steady state solution of 
Eq. ([6]) as a weighted ensemble of pure states in infinitely 
many ways, p ss = J2 k Pk\i>k)(ipk\, but only some of 
these are so-called physically realizable ensembles [13], 
i.e., there exists some continuous measurement of the 
environment such that the state is collapsed to \ipk) & 



proportion of the time pk- This has been referred to as 
the preferred ensemble fact. These ensembles are just the 
ones attainable through the quantum trajectory method, 
which gives a physical interpretation to the average gen- 
uine correlations studied above. 

Conclusions. — To define a measure of genuine multi- 
partite correlations suitable for the dynamical study of 
open quantum systems, we have exploited the fact that 
in many experimentally relevant situations the dynamics 
naturally give rise to a decomposition of the state into an 

0.8 




FIG. 3: Average genuine iV-partite correlations versus h/X for 
N = 100 and homodyne detection with two different detection 
angles. The angle 9 — 92.8° maximizes the correlations at 
h/\ = 0.5. 



ensemble of pure states. Even though the decomposition 
might not have any clear information theoretic interpre- 
tation, it has the advantage of arising from a well-defined 
physical process, as well as being well suited for com- 
putation through computer simulations. In a sense this 
measure generalizes the measure for average bipartite en- 
tanglement over quantum trajectories defined previously 
Q, and we have explored how going from bipartite to 
multipartite systems can give new exotic behavior in the 
dynamical evolution of average genuine correlations. Just 
as for the bipartite case we have seen that the measure 
can be sensitive with respect to the choice of unraveling, 
U, which leaves many questions open with respect to the 
significance of different detection schemes. Finally, we 
have shown numerical evidence that the proposed mea- 
sure of genuine correlations can serve as a witness for 
critical behavior in open quantum systems. 
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